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Abstract 

We consider a weakly interacting finite wire with short and long range interactions. The long 
range interactions enhance the Akp scattering and renormalize the wire to a strongly interacting 
limit. For large screening lengths, the renormalized charge stiffness Luttinger parameter -f^e//. 
decreases to Keff. < giving rise to a Wigner crystal at T = with an anomalous conductance 
at finite temperatures. For short screening lengths, the renormalized Luttinger parameter K^ff^ 
is restricted to ^ < Keff. < 1- As a result, at temperatures larger than the magnetic exchange 
energy we find an interacting metal which for K^jf. ~ ^ is equivalent to the Hubbard [/ — > oo 

2 

model, with the anomalous conductance G « ^ . 
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1. Introduction 



The anomalous conductance G ~ 0.7 x (2e^//i) discovered by Pepper et al. [l] and 
furtlier investigated by ,l6| is one of the major unexplained effects in quantum wires. 



Several theories have been proposed: Phenomenological theories 



polarization 



point contacts [1 
zigzag structures 



7|, l8|, Kondo effect |4|, spin 
formation of bound states [SH, unrestricted Hartree-Fock calculations for 

ferromagnetic 



Wigner crystal ferromagnetic spin coupling 



14j , the spin incoherent Luttinger liquid [15 



and the formation of 



a quasi-localized state ISj, however no consensus has been reached. 

For noninteracting spin unpolarized electrons, the conductance of narrow ballistic quan- 
tum wires connected to two (large) reservoirs is quantized in units of 2e^/h. An early 
suggestion was that the electron-electron interaction should modify the conductance for a 
Luttinger liquid 1^, 2o| as K{2e^ /h) where K is interaction-dependent. Using the method 
of Bosonization for weakly interacting fermions, it has been shown that by taking first the 
frequency limit — before the momentum limit q — » 0, the non-interacting leads modify 
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the metallic conductance to the limit C = ^ 

GaAs / AlGaAs in the lowest populated conduction band is a weakly interacting metal 
characterized by the Luttinger liquid charge K < 1 and spin Kg \ parameters. The 
presence of the unscreened, long range Coulomb interaction (in a one dimensional wire) 

we find that the Coulomb 



alters this picture. From our Renormalization Group study 
long range interaction enhances the weak Akp scattering channel and decreases the Luttinger 
charge parameter to -ft'e//. << 1 . As a function of the screening length, we can have either 
a strongly interacting metal (similar to the Hubbard U = oo model) or an insulating Wigner 
crystal. In order to capture both phases, we investigate a microscopic model of a wire of 
length L with a lattice constant a in the presence of a weak scattering potential. The lattice 
constant a is much smaller than the transverse width d of the wire, which controls the low 
energy excitations in the lowest band. The effective model in the lowest band is given by 
the renormalized lattice model. This is achieved by a Real Space Renormalization Group 
procedure, which replaces the discrete lattice model (lattice constant a) by the new lattice 
constant d = Integer x a, and renormalized interaction coupling constants. The effective 
model at the length scale d will be a function of the microscopic Fermi momentum kp (defined 
by the electronic density) and the effective umklapp momentum G*^^-* obeying the relation 
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G('^) . d = ■ a, where G^"^ is the microscopic umklapp vector. At the length scale d, the 
long range interaction is separated into two parts: the large momentum transfer, included in 
the effective short range Hubbard interaction, and the forward Coulomb interaction. The 



forward Coulomb interaction gives rise to Akp = G^^'sner) ^ _2n_ ogdllations where 
r^-e is the inter-particle distance. At T=0, a Wigner crystal ground state with a charge 
gap A is formed if the effective charge parameter obeys K^ff. < |. At finite temperatures, 

2 

comparable to the charge gap A, the conductance is given by G ~ ^. For short screening 
lengths, the interacting parameter Kejf. is restricted to | < K^ff, < 1 . As a result, we find 
that at finite temperatures which are larger than the magnetic exchange energy, the limit 
Keff, ~ I is equivalent to the Hubbard U = oo model. Therefore, at finite temperatures we 

2 

find the conductance is given by G ^ ^. 

The plan of this paper is as follow: In chapter 2, we present the interacting Fermion model. 



The renormalization effects 
mode formulation 
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br the finite wire will be investigated in chapter 3 using the zero 
27j. In chapter 4, we present the Fermion-Boson representation for 
the interacting wire. In chapter 5, we present the Renormalization Group (R.G.) analysis 
and show that the Renormalization effects of the effective charge interacting parameter 
Keff are controlled by the electronic density and screening length. In chapter 6, we use 
the renormalized interaction parameters to compute the effective zero mode Hamiltonian at 
finite temperatures for | < Keff. < 1. Chapter 7 is dedicated to the computation of the 
conductance at finite temperatures for | < K^jf, < 1. In chapter 8, we consider the case 
^ ~ Keff. and show that the model is equivalent to the incoherent Luttinger liquid which 
emerges at finite temperatures for the Hubbard f/ — oo model. In chapter 9, we consider 
the case Keff. < \ which at zero temperature gives rise to a Wigner crystal with a charge 
density wave gap A. In chapter 10, we present our numerical results using the experimental 
parameters given by |3i]. In chapter 11, we examine the effect of the Zeeman interaction. 
Chapter 12 is devoted to conclusions. Appendix A deals with the thermodynamics of the 
zero modes, and in Appendix B we present the calculation of the self energy for a wire of 
length L k, 10~^ at temperatures T f» 1 Kelvin. 

2. The model 

We consider an interacting wire at low electronic densities that has a finite width d. 
The geometric parameters in the quantum wire experiments are: the gate screening length 
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^ = 10~^m, the wire length L ^ 10~^m, the width d ~ the two dimensional carrier 
density ~ 2.5 x lO^^cm^^ and the electronic lattice spacing a ~ lO^^^m. Due to the 
width d, the single particle excitations are characterized by a set of electronic bands with the 
transverse quantization energies ^2 , r = 1,2, 3.... The gate voltage is such that, at the 
temperatures considered in the experiment, only the lowest band is populated. In order to 
describe the low energy physics, we project the microscopic Hamiltonian iJ" (the microscopic 
model defined at the lattice scale a) into the lowest band. As a result, the effective one 
dimensional Hamiltonian H with the energy cut-off characterized by the transversal energy 
separation momentum cut-off A = ^ preserve the original form of the 

microscopic Hamiltonian if": 

+U T 'ip'i {nd)ip i{nd)ip^ {nd)ip'^ (nd) 

n 

+ EE E T ^Pt{nd)Mnd)V^''\\nd-n'd\)tP+{n'd)tP^,{n'd) (1) 

where ~ 2m*d'^ effective hopping at the length scale d, U is the projected repulsive 

Hubbard interaction, which also contains the effect of the Coulomb interaction obtained 
by projecting out states with a lattice spacing in the interval a — d. V^'^\\nd — n'd\) = 

, ^ , ^ = is the effective long range Coulomb interaction defined for 

distances x > d and ^ is the screening length. The Hubbard model is characterized by the 
particles - holes charge and spin excitations: K < 1 (charge) and > 1 { spin ). We 
consider the situation where the Fermi momentum kplVc) satisfies the condition AkpiVc] 
= Te-e > d, suggesting that the umklapp interaction is negligible. According to 
the density is expressed in terms of the external gate voltage Vq: kpiyc) = ^neiVc) = 
~"(Vg — V^'^), where V^'^ is the gate voltage at which the wire is pinched off and ^ is the 
effective capacitance. The results reported in ]^ show that the conductance decreases with 
the lowering of the gate voltage, suggesting that the umklapp interaction is significantly 
enhanced. The low energy properties of the model will be investigated using a combined 
method of Bosonization and R.G. theory. At finite temperatures, the exact description of 
the electron excitations requires the inclusion of the zero — modes operators. 

3. The representation of the electron operator for a wire of length L 

The electron is represented as a product of two operators, a Bosonic one (this is the 
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standard Bosonic representation for spin-charge excitations) and a Fermionic-zero mode 
operator, which carries the electron number (electrons with spin up or spin down that are 
added or removed from the Fermi surface). The electron operator ip'^{x) is restricted by 
the momentum with a momentum cut-off [A, —A] around the Fermi surface and is given in 
terms of the right Ra{x) and left L„{x) components: ipai^) = e'^''^^ R„{x) +e~^''^^L„{x) with 
the Fermi momentum kp = /cf(Vg). We replace the right (left) mover fermion by a product 
of a fermion operator Fji^^{x) {Fl^^{x)) and the boson one e'-^^^R.^ix) ^g«v^0L,<T{a;')); 

R^{^jc) = P^^ic^R,.^ii2n/L)(Nn,^-i/2)x^iVA^enA-) = F^^^(x)e'^®«'-(") (2) 
V 27r 

L^(^x) = \ ^ ^-iaL.cr ^-iCi-n / L)(NL.a-ll'2)x ^-iVA^eL.cr(x) ^ Fl a{x)e^'^^^^'^'''^^^ (3) 



where e^^^'^R^-^i'-^) and e~*^^®^.<'(^) are the standard Bosonization formulas used in the liter- 
ature. FR fj{x) and Fl^„{x) are the zero mode fermion operators that can change the number 
of particles and are crucial for enforcing the Fermi statistics. This operators are defined with 
respect to the non - interacting ground state \F > (the Fermi see characterized by the Fermi 
momentum) . The electronic Hilbert space excitations above the Fermi see are given by the 



states 
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26l |: rrig > ^iN^y, m'^ > where > , > are integers which specify 

the number of Bosonic quanta (particles -holes excitations) with a momentum q = > 0. 
Nr^o- , Nl^„ represent the change of the total number of electrons in the right and left ground 
states. The formal proof that relates the electron operator to the zero mode operators is 



given by the Jacoby identity 
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26| are given 



27l |. The Bosonic - particle hole excitations 
by: 0R,o-(2;) = 0cr(a^) — ^a{x), 0L,o-(x) = 6o-(x) + ^a{x). The zero mode Fermion excita- 
tions are given in terms of the zero mode coordinates aFi.^a,OiL,a and their canonical conjugate 
fermion number operators iV^ o- , N^^^ ( cr =1, |). The physics of the zero modes is described 
in terms of the charge operator Qc = J2a=i,^WR,cT + ^L,a] and the magnetization operator 
Qs = [(A''/jo-=t + ^L,o-=t) " {^R,u=i + ^L,a=i]- The canonical conjugate variables to the 
charge and magnetization are given by the charge coordinate d = J2a=i,^[(^R,cT + CiL,a] 
and the magnetization coordinate as = [(«R,o-=t + «l,(t=t) ~ («i?,o-=j. + aL,<T=j,)]- The 
zero modes obey the commutation rules : [aR,a, Njiy] = iS^y, [—aL,a, Ni^a'] = i^ay and 
[an^a, Niy] = [oiL,a, Nr^] = 0. 

4. The model Hamiltonian in the Boson-Fermion representation 
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The Bethe Ansatz formulation [28| and eqs. (2 — 3) allows us to map eq. (1) into a charge 
and spin interacting model. The mapping is a function of the Hubbard interaction strength 
U = J- and the electronic density Ue- The Hamiltonian is controlled by the charge parameter 
K = K{U,ne) , spin parameter Ks = K{U,ne), umklapp interaction g = g{U,ne) = gA."^, 
spin backward scattering parameter gs = g{U,ne) = (jsA"^, Fermi velocity vp, charge density 
wave velocity v = v{U,ne) and the spin density wave velocity Vs = Vs{U,ne). The Hamil- 
tonian in eq.(l) is replaced by: H = H^^^ + + H^^^^\ The first two Hamiltonians 
H^^'^ + H^^^ represent the particle hole excitations and if("=*') represents the zero modes. 
The charge excitations if"^*'(0, $; d, Qc) are given in terms of the Bosonic fields 9 = ^^J^-^ , 
$ = "^^J^^ and the zero mode Fermionic fields a,Qc- 

rL/2 X 1 

J-L/2 2 2K 

2% ^ 

-g cos[V8nQ{x) + a + (4A;j.(Vg) + -rQc)^]]] 

ij 

p2 .L/2 fL/2 p2 p2 

+ / / dxdx'd^Qix)\ , - - , ]d.,'Q(x) 

TlKo J-L/2 J-L/2 ^^x-x'f + d? ^ {x - x'f + ^ 

(4) 

where v is the charge velocity vK = vp = and K = K(U, Ue) is the charge inter- 

action parameter. The last term in equation (4) represents the forward part of the long 
range interaction given in eq.(l) with the screening length ^. The long range interactions is 
controlled by the coupling constant 7 = ^ ■ where c is the light velocity and kq = 13.18 
is the dielectric constant for GaAs. The strength of the umklapp interaction g is determined 
by two parts: the short range Hubbard U repulsive interaction and the large momentum 
transfer of the Coulomb interaction obtained after the projection. The spin density wave 
excitations are given by the Hamiltonian H^^^{Qs,^s]C(s,Qs) with the spin density wave 
operators: 6^ = = 

/•L/2 X 1 
H:^\Qs, as, Qs) = Vshl / dx[^{d,<^s{x)y + :rTr{d.Qs{x)f 

J-L/2 I IKs 

+gsCos{V87!-Qs{x) + + —Qsx)]] (5) 

where Vs = < v is the spin wave velocity, > 1 is the spin stiffness and gs is the 
spin density wave coupling constant. Next we present the zero mode Hamiltonian if("-=o) 
using the normal order notation: : if("=°) := ff("=o)— < F\H^^^^^\F > where \F > is 



the unperturbed Fermi surface at zero temperature given in terms of shifted operators, 
NR,a = Nn,a+ < F\Nr^„\F >; NL,a = NL,a+ < F\NL,a\F > (see Appendix A). The normal 
order, zero mode Hamiltonian takes the form: : if :=: Hq"'~^^ : + : H^^t^'^^ '■■ 



: Hf' ■■= ^ [Nl,.^, + Nl,., + Nl.^, + iV^^J (6) 



: H^:r'^ := n(^)(L)[(iV«,.=t + Nl,.=^) + {Nr,^=i + NL,a=i)? 

The zero mode couphng constants obtained from eq.(l) are given by the renormahzed charge 
backward interaction u^^\L) = ^ff- C~[^^ ) and the backward spin interaction u'^^^L) = 
^^ C'k^" )• temperature and L ^ oo, Kg flows to 1 and the backward interaction 

vanishes. The function F(^^) = /o^[^^^S±ll - /o^fV^^^i jg ^j^g Fourier 
transform of the long range screened potential. At flnite temperatures, the Fermi energy is 
shifted by 6fio{T) which modifles the zero mode Hamiltonian: 5H^'^""^^ = 6fio{T)[{NR^„=i + 

5. The Renormalization Group equations 



One of us 



30(1 has developed an R.G. method which is applicable for the Hamiltonian 



representation. This method has been used [2^ to derive the R.G. equations for the unbiased 



Sine-Gordon in the presence of long range interactions controlled by the coupling constant 

=2 1 
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he kq ' 



1-1/2 X 1 , 

J-L/2 I IK 

+ / / dxdx'dMx)\ , =\d^,Q{x') 

TTKo J-L/2 J-L/2 J(^^^'\2j_^2 



TlKo J-L/2 J-L/2 ^(^x-x'y + d^ 

(8) 

In the absence of the Coulomb interaction the model is equivalent to the classical two 



dimensional Sine-Gordon model. According to [33|, |3J] the model is gaped for K < 
n = 1,2.... The long-range interaction modifles the results and drives the model to a gaped 
phase for any value of K. We have extended the R.G. calculations for the biased Sine Gordon 
model 5fcos[\/87r9(a;) + a + {Akpiyc) + ^Qc)x] given in eq.(4). We obtain the new R.G. 



equations as function of the bias and the screening length ^ for the differential momentum 
shell dl = -!f. 

{AkFiVc) + —Qc) ^ i^kFiVc) + —Qc)e' (9) 
dl '"^^^ ^/(l + 7(3n)MH(/)) 4(l + 7(-^)M,(/))^ ^^"^ 



dl 8(l + 7^M«(/)) 



(11) 



(12) 



dl 4(1 + 7^M«(0) 

where Mr{1) is the difference of two Bessel functions Kq{x): Mr{1) = 2(Ko[e^'] — ii'o[| -e"']) 
The solution of the R.G. equations depends on the initial values of the interaction parameters 
gnil = 0) , Kji{l = 0) and the ratio ^. We will study the case where AkpiVG) < ^• 
In order to compute the scaling functions, we need to determine the relation between the 
logarithmic scale I and the voltage Vq- Based on the experimental observation [sl we have 
a perfect conductance for a particular gate voltage Vq'^ for which the umklapp interaction 
is negligible. This will happen if AkpiVQ^) corresponds to the momentum ^. For this case 



we find an oscillating behavior for the Sine Gordon term: g cos[v 87rG(a;) + d + (4fci?(V^ ) + 
^Qc)x] = 1)" cos[-\/87r6(x) + a + {^Qc)x] and can ignore the umklapp contribution. 
For lower gate voltages Vq < Vq^ the situation is different. Following 31| we do not 
neglect the umklapp interaction for Vq < Vq\ instead we compute the effective coupling 
constant at the length scale / = /(Vg). This length scale / = liVc) is determined by 
the equation 4fc^(V^G)e'^^'=) = ^kpiV^^^) and is given by 1{Vg) = At this 

length scale, the renormalized umklapp interaction alternates in sign g{—l)"' and can be 
neglected if the Sine-Gordon coupling constant is small. Using this procedure we substitute 
the function 1{Vg) into the R.G. equations and find the renormalized Luttinger parameter 
as a function of the gate voltage Vq- Since the wire has a finite length L we stop the scaling 
when we reach, the value / = minimum[l(yG),lL] where II = log[^]. In the presence of 
the Coulomb interactions the Luttinger parameter Kfi{l{VG)) is replaced by the effective 
parameter KeffXliyc)), computed from the R.G. equations (9 — 12): 
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The effective interaction parameter Kef/XK^c)) decreases monotonically with the decrease 
in the density and exhibit a maximum for densities where KeffXK^c) ~ |- The charge 
density velocity is enhanced to f = ^ f^{i(VG)) ' ^^^^ screening ratio approaches ^ = 1, 
the Coulomb renormalization is absent and KeffX^iVG)) = Kr(1{Vg)). In figure 1 we have 
plotted -j^ — Ii{Vg)) ^ function of the gate voltage. Following [3i], we have used for the gate 
voltage Vq^ the value Vq^ = —5.1 volt. We observe that — {i{Vg)) ^ minimum for volt- 
ages that corresponds to the region where the 0.7 feature is seen. Since the compressibility k 
is proportional to the square of the Luttinger parameter k oc ii'gjj (/(Vg)), we conclude that 
a maximum in the compressibility suggests the formation of a gap. (Since the compressibil- 
ity is proportional to the derivative of the renormalized chemical potential /iR(VG', /(Vg)), 
^ = {ne(yG)ydn^(VG)[f^R(yGj{VG))] wc cxpect also a minimum for the derivative.) 

6. The effective Hamiltonian i < fs:e//.(/(Vb)) < 1 for L > ^ 

Using the dependence of the Fermi momentum kp{VG) < kpiVQ^) on the gate volt- 
age Vg < Vq \ we find that the umklapp interaction and the Luttinger parameter are 
renormalized. Following the analysis from chapter 5, we find that at the length scale 
K^g) = ^'^di '^ikp^ya) ] renormalized umklapp interaction is negligible (/(/(Vg)) ~ and 
the renormalized velocity is ^ ff{i{VG) ' ^'^^^ L > C,, the effective Luttinger parameter is 
restricted to | < K^f/XK^G) < 1- The renormalized Bosonic Hamiltonian is given by: 



H / dx dx'dx&Rjx) \—^ =]dx'QR(x') 



(14) 



H:^\e,,^,;a,,Q,,l{VG)) = Ki?(^(^G))[/ , dx[ ^^\ g.1>.,^(x))^ 

+^7Iw77T7t(^-®^'«(^))'] + 9s{1{Vg) cos(v^e,,^(x) + + ^Q.xe'^^^))]] (15) 

Since Ks{1{Vg)) > 1, the Sine-Gordon scaling shows that gsiK^G)) is an irrelevant coupling 
constant which decreases with the increase of 1{Vg)- The renormalized zero mode Hamilto- 
nian will depend on the renormalized coupling constants given by the R.G. eqs. (9 — 13): 

: H^''='\1{Vg)) :=: Ht'\liVG)) : +Htf\l{VG)) : (16) 



The first term : H^^ °\1(Vg)) : represents the non-interacting part: 

: Ht'HliVa)) := ^[^l.^T + ^L^T + ^l^-i + ^l.^J (17) 

The second term represents the interactions ■.Hj;^^'^\l{VG)) '■, given as a function of the 
charge operator = [(A^ii,a=T + ^L,a=T) + (^ii,a=i + ^L,a=i)] = 4+ < F\Q^\F > and the 
magnetization operator Qg = [(-/Vr,(t=t + NL,a='[) - {NR,a=i + NL,a=i)] = Qs+ < F\Qs\F >. 

■ Ht:'\l{VG)) := t^MVgWc - Vs{1{Vg))QI (18) 

where vMVg)) = ^iC-^^^) + 7{^)Fi^J^, iJ^)] the renormahzed backward 
charge and magnetic interactions ?7s(^(Vg)) = ^f£- ( ^ k^"'1i(v^)^ ) ■ Both terms are a function 
of the screened Coulomb interaction F( ^uva) i de'(^G) ) gi"^®^ by: 



L c y[i+(^)^]+i . ■ V[i+(^)^]+i . 

At finite temperatures the effect of the e-e interactions replaces the non-interacting ground 
state \F > with a shifted Fermi surface given by the renormalizcd ground state |G > . In 
Appendix B we find that the single particles states e(n) are shifted up in energy by the self 
energy ST,{Vg,1{Vg))- In Appendix B we have computed the self energy ST,{Vg,1{Vg)) at 
low temperatures T which are higher than the spin exchange energy, KbT > r]s{l{VG)) < 
G\Q^g\G >= KbT* . > is the renormahzed Fermi Surface which replaces the non - 
interacting Fermi surface \F > and T* is a temperature of the order of 0.05 Kelvin. For 
temperatures T > T* the self energy is given by 6J:{Vg,1{Vg)) ^ 2r]^{l{VG)) < G\Qc\G > 
(see Appendix B). The effective zero mode Hamiltonian is replaced by: 

+5E{Vg, l{VG))[NR,a=^ + NR,a=l + ^L,.=T + ^L,.=J (20) 

7. The current for the interacting region | < Xe//.(^(Vb) <l,T>T* 

For finite values of 1{Vg) the spin density wave coupling constant gs{l{VG) and the spin 
density wave velocity Vs{1{Vg) — luva) k {UVg) ~ "^{K^g) are both small. At temper- 
atures T > T*, we replace the interacting zero mode Hamiltonian with the effective zero 
mode Hamiltonian controlled by the self energy ST,{Vg,1{Vg)) given in eq.(20). In order 
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to compute the current, we include the reservoir Hamiltonian Hues controlled by the drain 

(0) _ (0) 

source voltage V = ^RO/ht . 

° e 

Hb.s = ^ E mL,a - NR,a)] (21) 

The partition functions in the presence of the reservoir is given by: Z — 

T^[e-P--H^''=°Hl{VG))efr-e-P--Hne.--j = 2;[g-/3:i/("=°):g-/3://|/^{ :]_ rj.^^ gg^f ^^^^^^ ^^^^g ^g rCplaCC 

the reservoir Hamiltonian i^ijes by an effective reservoir : 

H'£- = Hues + E [^^(VG: l{VG)){NL,a + Nn,a)] (22) 

The static conductivity is computed using the non-interacting zero mode Hamiltonian ifo"~°^ 
given in eq.(17) and effective reservoir H^' given by eq.(22). The current is obtained from 
the derivative of the zero mode coordinate a (see chapter HI ), I — Using the 

Heisenberg equation of motion we obtain the current operator. 

o'=T,J. 

The thermal expectation function is obtained with the help of the partition function Z. 

I = r,[e-^<"'e-'^<^/][Z]-i (24) 
Following Appendix A we obtain: 

- \[JF.D.[ Jl—=, J 

-/F.o.l -j^^ J J l^^j 

where €1(171) and e/j(m) are the single particle energies and 2ni?(VG) is the discrete bandwidth 
introduced in Appendix A. We include a small single particle broadening which will allow 
us to replace the discrete sum eR^L{m) by a continuum integration variable e. Performing 
the integration with respect the energy variable e and expanding with respect the voltage V 
gives the conductance G — ^: 

n 2e , 6 + ^S(l^G,/(KG)) + f -^MT) ^ ^ 6 + ^S(K;, /(\/c)) - f - ^/Xo(r) ^^ 

^ " wi-eH.o)^'^^""f ^"^""f 

2e%„ -6p(Va) + ^S(Vb,/(VG))-^/xo(r) ^ ^ 6p(Vg) + ^S(Vc, /(Vb)) - ^/xo(T) ,, 
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We observe that the self energy determines the conductance through an effective chemical 
potential. The bottom of the bandwidth —eFiVc) is replaced by —eFiVc) + 5S(Vg, ^(Vg)). 
This allows to introduce the renormalized effective chemical potential /x_r(Vg', /(Vg)) = 
eF{VG)-6J:{VG,l{VG)). 

8. The strongly interacting region KeffXK^c)) ^ ^ , T > T* - The effective 
U = oo Hubbard model 

When KeffXK^c)) ~ | ^-nd T > T*, one obtains an incoherent Luttinger liquid which can 
be mapped to the Hubbard model U ^ oo. (When f/ — oo the interaction Luttinger pa- 
rameter is given hj K ^ and the spin excitations which are of the order ^ can be ignored.) 
This limit U ^ oo has been considered in the past j29]. In this limit the following con- 
straints must be obeyed: il>t-i(x)il>„=^(x)+'il't-i(x)il>„= \ (x) = 0, 1. Using the constraints, we 
have found the following representation 29| for the electron operators: ipa{x) = 6o-(x)^(x), 
ipti^) ~ '^~^{x)b^{x) where is the electron charge operator and 6o-(x) are the hard core 
boson for the spin excitations. They obey the constraints: 6^^^ (x)6o-=| (x) +fe^=i {^)b(T=l = 
\t'+(x)\l'(x). In one dimension, this model has been represented in terms of the Bosonic elec- 



tron operators 



291] 0e and $e and Spinon operators B^, and The constraint is imposed on 
the electron density: pe{x) = p^=t(x) +p^=|(x) = -^[d^Q^=^{x) + d^Q^=i{x) = ^9:r6e(x). 
The canonical conjugate momentum is given by: dx^e{x) = ^[dx^a=i{x) + dx^a=i{x)]. For 
non - interacting electrons we have the commutation rule [Be(a;), 9^.$e(y)] = ih5{x — y). 
Due to the exclusion of double occupancy, the electronic density is reduced by a factor of 
two (in comparison with non - interacting electrons) and the commutator is modified to: 

[6e(a;), dx^e{y)]constraint ~ ^M{x - y). 

The Hamiltonian for the f/ — > oo case (away from half filling ) is given in terms of the 
fields ee(x), $e(a;): H, = ! dxvh[{dx^,{x)f + \{dxQe{x)f]. 

At finite temperatures, the spinon Hamiltonian i/, = / dx^[{dx^s{x)f + {dxQs{x)f] ~ 
is negligible ( e,(x) = ^[e<,=t(x) - Qa=i{x)] and $,(x) = ^[<l>^=|(x) - <l>^=^(x)]). If we 
inject an electron with a given spin at one lead, we will detect on the other lead a charge 
with an arbitrary spin. The effect of voltage difference V between the leads, is included 
into the calculation through the reservoir Hamiltonian f dx[dx^a='\{x) + dx^a=i{x)] = 
^ / dxdx^e{x). Using the Heisenberg equations of motion with the modified commuta- 
tor we obtain the electronic current operator [37| Je = ^^[dx^a=^{x) + dx^a=i{x)] = 
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-^dx^e{x). The extra factor of | which appears in the current operator Jg is due to the 



commutator 



.) J Constraint 



37|. Therefore, the conductance is reduced to G 



9. The conductance in the Wigner crystal limit KeffXK^c)) < \ 

For large screening lengths the effective Luttinger charge stiffness i^e//.(^(^)) decreases 
below KeffXK^G)) < I at low temperatures. Under these conditions, the R.G. analysis 
reveals that the alternating umklapp coupling constant g{—l)'^ generates a gap at T = . To 
investigate this region, we introduce two sub - lattices for the even and odd sites. We replace 
the Bosonic fields by the even and odd combinations : 0_ {y = 2nd) = Q(^-2"°')~^^-(2"+^)'^) 
and 6+(?/ = 2nd) = ^ integrate out the antisymmetric field Q-{y = 

2nd) = ^ Q^Yid obtain an effective Hamiltonian for the symmetric field 

©+(?/ = 2nd) = _ 'jj-^g effective Hamiltonian has a set of new coupling 

constants \g^f}^\ ~ ^ < (sin[V47r0-(y)])^ >. 

i7,"^°(e+, $+) ^hv{l{Va))[j"] dy f^f^-l^^''^\ dy^^{y))' + ] .A dyQ^v))' 
+9'^l,RmG)) cos(v^0+(y))]] (27) 
where the new couphng constant Qnew rQ'^g)) obeys the R.G. equation. 

rf^ieU(0_.,(2)_^(^^j^_2^^^^^(^)j (28) 



dl 

^(2) 



inew, 



This equation shows that g\^^^ ^{l) is a relevant coupling constant for KeffXK^c) < |. As a 

1 

result, a charge gap A ^ A{g';^JJ^^'"^''ff-'-"-'^G))-^) ^ill open. When KeffiK^c) < |, we obtain 



from eq.(28) that at T = the expectation value of the phase Vl6vr < Q+{x) >= tt will give 
rise to a Wigner crystal order p(x) ~ constant+cos[4:kpx+-^]e~^^^®^^^~^^^^''-^'''^^ . (Expand- 
ing the cosine term in eq.(27) around the ground state < <d^{x) >^ shows that the charge 
density wave has a gap A. This gap suppresses the fluctuations e~5<(®(^)~<®(^)>)^> ^ 
and stabilizes the Wigner Crystal order at T = 0.) In order to evaluate the effect of the 



charge gap A on the electronic spectrum we map [35|, |36| the Bosonic charge Hamilto- 
nian to a spinless Fermion for KeffXK^c) ~ |- We introduce a two component spinless 
Fermion: x~'"(x) = [xii^) ^ X2{x)]'^ = yA[e*v^e+(a;)^ g-^^/4;F0+(x)j + _ ^ result we find for 



KeffXK^c) ~ I fhat the Hamiltonian in eq. (27) is mapped to a spinless Fermion model: 

H,,p = J dx[hv{liVG))xHx){-tdxas)xix) + ^'''''"'''f^''^\ x^{x)a,x{x))'] 
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^ J dx[hv{l{VG))x^ix){-id^a3)xix) + g^^^^ j^iliVc)) < X+(x)(Tix(a;) > X+(a;)(Tix(x)] 

(29) 

where as and ai are the Pauh matrices. As a result, we have a gap 2 A between the lower 
band and the upper band given by the self consistent solution: A ^ A = ^ieL _r(^(^g)) < 
X~^{x)o'ix{x) >■ The energy difference A between the Fermi energy and the top of the 
lower electronic band will affect the conductance through the Fermi - Dirac function. For 
this case, the self energy is replaced by the gap A and the conductance is approximated by 
G ^ ^[1 — fF.D.ii^)]^ for ^bT > A conductances is given by G ^ ^. 

10. Numerical results 

We have used the experimental relation between the Fermi momentum and the gate 
vohage Vg given by Kf{Vg) = frie(^G) = ^{Vg - V^'') , V^^ ^ -5.52volt, ^ = 1.2 • 
10^{Volt ■ meter)^^ [3] to compute the conductance in figures 2 and 3. In figure 2 we 
have considered a typical screening ratio | = 10 and plotted the conductance for a varying 
range of temperatures 1 — 3 Kelvin. Figure 3 shows the conductance at a fixed temperature 
T = 1 Kelvin for a different screening lengths. We observe that for ^ = 1, the Coulomb 
interaction is completely screened and the 0.7 feature is absent. In figure 4, we plot the 
dependence of the self energy 5S(Vg, / (Vg)) on the gate voltage Vg- We observe that at low 
densities the free energy has an extremum at a finite density. The renormalized chemical 
potential yUR(K;, /(^g)) = eF(Vb) - 5S(Vb,/(Vb)) vanishes at the voltage > V^^. For 
V*^ <Vg <Vq, the renormalized chemical potential is negative, indicating the formation of 
a charge density wave gap at T = for Kf,ff {l{yG) < \- The derivative of the conductance 
and chemical potential are related to the inverse compressibility: '^^^^^ oc "'^-^(^^^(^g)) (x 
o.,, \ The 0.7 anomaly is translated into a minimum around Vg = —5.49 volt for 

the conductance derivative ^'^^^-^ and the compressibility k(Vg) which is proportional to 
the inverse square of the effective interaction parameter ii'e//. (/(Vg)) shown in figure (1). 
Therefore, we have the confirmation for the formation of a charge density wave gap for 
Keff.{l{yG))i\ at zero temperature. In figure 5, we plot the function '^t^R^'^^^^^^)) ^ ^YhSs 
function has a minimum at the voltage Vg = —5.49 volts, which corresponds to the 0.7^ 
structure observed for the conductance graph. 
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11. The effect of the Zeeman magnetic field 

The Zeeman magnetic field introduce a bias term {2{kp(yG) — {kp(yG))x into the last 
term in eq.(5). Expressing the bias in terms of the magnetic field B\\ we find, 2{kp(yG) — 
ikUVc) = 2kp{VG){^l + ^-^l^) ^ ^kp{VG){i^) ^^^^^ = 
energy and = g\\fiBB\\ is the Zeeman energy. As a result the spin part Sine-Gordon term 
vanishes since Ks{l) > 1. As a result the spin wave velocity = ^ is further reduced. 
For large magnetic fields > 1, the wire will be polarized and we will have only one 
propagating channel with the conductance G ~ 0.5 x (2e^ /h). 

The effect on the charge density wave Hamiltonian will be to replace AkpiVG) in eq. (4) 
by: 2[kl{VG) + 4(Vb)] = ^kF{VG)[^l + ^ + - ^]\. For |j < 1 we show that 
the perfect conductance in the absence of the Zeeman magnetic field computed at the gate 
voltage Vq is shifted to a larger gate voltage y^-^'^'^^'^^ in the presence of the Zeeman 
field: 2[4(^g"^'~) + k^piy^-^™"")] ^ AkF{y^-^™'')[l - \{^?] = ^kpiVS) = f • 
This formula shows the shift in the perfect conductance from 4/ci?(Vg) = ^ to a larger gate 
voltage y^o-Zeeman) ^ yQ gj^^^ /cp ( Vq""^*^^™""^ ) ~ T—fr^T^- This rcsult is in agreement 
with the experimental observations [2j. 

The conductance at finite temperatures T > T* will be given by replacing the self en- 
ergy in eq. (26) with a new self energy computed in the presence of the magnetic field, 
ep,^=^{VG) = cfIVg) + ^ and eF^fj=i{VG) = ep{VG) — The results for the conductance 
are shown in figure 6. We show three graphs: the first graph (thin line) represents the 
conductance in the absence of the magnetic field and the other two graphs represent the 
conductance for the magnetic fields 5 = 3 Tesla and i? = 10 Tesla. We observe the shift 
of the conductance to higher voltages with the increase of the magnetic field . 

12. Conclusion 

We have presented a model which explains the conductance anomaly at finite temper- 
atures as a function of the gate voltage. Due to the Coulomb long range interactions a 
weakly interacting electronic system can fiow to the strong coupling limit Keff.(l{^G)) < \- 
When the screening length is not too large, the Luttinger stiffness is restricted to | < 
Kf.ff {HVg)) < 1 . As a result, the conductance of an infinite wire is perfect at zero tem- 
perature. At temperatures larger than the magnetic exchange energy T > T*, we have an 
incoherent Luttinger model. For Keff.{t{yG)) ~ \ the interacting wire is equivalent to the 



15 



2 

Hubbard U oo model with the anomalous conductance G ^ ^ . 

h 

For large screening lengths the interacting charge stiffness decreases to ii'e//. (^(Vg)) < \ 
. As a result we find that at zero temperature we have a Wigner crystal with a charge 
gap A. At finite temperatures the formation of charge density wave gap gives rise to the 
anomalous conductance G ~ ^[1 — /F.D.(i^T)]- Following [2] we have investigated the 
effect of the magnetic field. We have shown that the magnetic field shifts the region of the 
perfect conductance to higher voltages. 



Some of the concepts used in our work are common to other theories 38, 39 1 . The long 

n n J 6 

range interactions have been introduced by |12l . 1141 ]: in the present paper we show that, 
by varying the gate voltage and screening length we obtain either a strongly interacting 
metal or a Wigner crystal. Other theories use a weak scattering potential [81] or charge 
localization j^,|4^ and are consistent with our picture. In our view the origin of the scattering 
potential (microscopic or phenomenological) is not crucial! The crucial effect is that any 
weak scattering is strongly enhanced by the long range interactions! Our findings show 
that, due to the long range interaction, any negligible scattering potential is enhanced and 
eventually can drive the system to an insulating regime. It is the interplay between the 
screening length, gate voltage and temperature which gives rise to the conductance anomaly. 
One of the popular theories is based on the the Kondo model {4, 9]. The Kondo picture 
dictates that the anomaly should be observed above the Kondo temperature. When the 
temperature is lowered below the Kondo temperature, the conductance is restored to the 
universal value. This picture is consistent with our theory in the following way: If the strong 
coupling regime KeffXK^G)) ~ | is reached in the metallic phase, we can use the Hubbard 
U ^ 00 limit, which is the basis for deriving the Kondo model. The Kondo physics emerges 
for finite exchange coupling J oc ^ . In our case, the anomalous conductance is observed at 
finite temperatures T > T* which is comparable to ^ where the Kondo picture emerges. 
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The authors acknowledge the financial support from the CUNY Collaborative Grant award 
for the year 2007-2008. 



16 



0.016 r 

■ 

O.OR 

■ 

0.012. 

ir 

0.010 - 

■ 




-5.50 -5.45 -5.40 -5.35 -5.30 -5.25 -5.20 

FIG. 1: The effective parameter — \, , that is proportional to the inverse compressibihty is 



plotted as a function of the gate voltage / = 1{Vg) for, L = 10^^ meter and the screening ratio 
1 = 10 



G 




V_{g} 



-5.5 -5.4 -5.3 -5.2 -5.1 

FIG. 2: The Conductance G in units of ^ as a function of the bias gate voltage / = 1{Vg) 
for the temperatures T = 1. Kelvin (upper hne), T = 1.25 Kelvin^ T = 1.5 Kelvin,T = 1.75 
Kelvin,T = 2.0 Kelvin,T = 2.25Kelvin,T = 2.5Kelvin and T = 3. Kelvin (the lowest line) for 
umklapp parameter g{l = 0) = 0.05, K(l = 0) ^ 0.98 ,L = 10^^ m and screening ratio ^ = 10 
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FIG. 3: The Conductance G in units of ^ for four screening ratios | = l(upper line), | = 1.1, 
i = 1.3, i = 1.5, i = 2., ^ = 3., ^ = 5., ^ = lo,i = 50. and ^ = 100 at temperature T = 1. 
Kelvin length L = 10~^ m for the interactions parameters gR{l = 0) = 0.05, K{1 = 0) ~ 0.98 
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FIG. 4: The shift in the chemical potential 6T,{Vg, liVc)) for screening ratio | = 10 at temperature 



T = I.Kelvin length L = 10 for the interactions parameters gii{l = 0) = 0.05, K{1 = 0) ~ 0.98 
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d)U/dVo 




-5.50 -5.45 -5.40 -5.35 -5.30 -5.25 -5.20 

FIG. 5: The derivative of the chemical potential '^^^^-^y^^^'^-'-* for screening ratio | = 10 at tem- 
perature T = I.Kelvin length L = 10~^m for the interactions parameters gR{l = 0) = 0.05, 
K{1 = 0) « 0.98 
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FIG. 6: The effect of tlie magnetic field on the conductance in units 2e^/i. The first graph 
represents the conductance for zero magnetic field, the second graph represents the conductance 
for a magnetic field B = 3 Tesla and the third graph represents the conductance for the magnetic 
field i? = 10 Tesla. The other parameters were: screening ratio | = 10, temperature T = 1 
Kelvin, length L = 10"^ m, d = 10"* m, gR{l = 0) = 0.05 and K{1 = 0) « 0.98 



22 



Appendix-A 

The non-interacting Fermi surface at T = is given by the state |-F >, which is con- 
structed from the vacuum |0 >: \F >= n.=Ta[n-"/Ht) ^^(^' ^) n^^Cf^ ^^(^' >" 
We introduce the notation Nji^„ A^l.ct for the normal order at zero temperature: 

NR,a= E R^{rn,a)R{n,a)- ^ < F\R+{m,a)R{m,a)\F >= Nr,^- < F\Nr,^\F > 

(30) 

-npiVc) -npiVa) 

E L+{m,a)L{m,a)- E < F\L+ {m,a)L{m,a)\F >^ N^^,- < F\NL,a\F > 

np{VG) np{VG) 

(31) 

The presence of a reservoir with two chemical potentials ^r and //l is described by the 
reservoir Hamiltonian: 

npiVa) np{VG) 

Hrcs^/J'R E R'^{rn,a)R{m,a) + E L'^{rn,a)L{m,a) (32) 

-npiVa) -np{VG) 

At finite temperatures, the Fermi surface is shifted by SijLo{T) and is given (for the one 
dimensional case) by: Sfj,o{T) — eF(VG)f^(7^f^)^)- The temperature and the reservoir 
modifies the number of fermion in the thermal ground state to < Nr^^{Vg, Hr,T) > and 
< NL,a{VG,I^L,T) > given by: 

< J^rA^g,IJ'R,J-)>- 2^ If.dX 7— J {66) 

np{VG) 

< NL,a{VG,^L,T) >^ 2^ fp.DX 7^ J (34) 

np{VG) 

The expectation value of the normal order operators will be given by: 

< NLAVG.^J^L,T) >=< NlAVg,^^l,T) >-< NL,^{Sfio{T),T) >; 

< NR^^{VG,tiR,T) >=< NR^^{VG,fiR,T) >-< NR^^{6fio{T),T) >. 

The effect of the self energy 5T^{Vg, 1{Vg),T) will be taken in consideration by substituting 
in the previous equations : /ir ^ /ir — (5E(Vg, 1{Vg),T) and i^l ^ /J'L — S'EiVG, 1{Vg),T) . 
Appendix-B 

The purpose of this Appendix is to compute the self energy 5T,{Vg,1{Vg),T) for the 
following model: 

: H^::'\i) ■■= vmvg))qi - vsmG))Qi m 
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where Qc — (5c+ < F\Qc\F > is the charge operator and Qg — Qs+ < F\Qs\F > is 
the magnetization operator. We observe that the zero mode component of the Hamilto- 
nian commutes: [Hq"~^\1) , Hj;^^^\l)] = 0. Therefore, at finite temperatures, the parti- 
tion function = Tr[e~^'^*" °^ '•'-'] can be computed exactly. Our goal is to compute 
the charge current / = e^, which is given by the commutator [a, H^"'^^'>[1)]. We will 
limited ourselves to finite temperatures such that the exchange energy is smaller than 
the thermal energy and therefore, can be ignored (for long wires the spin stiffness ap- 
proaches Kg 1 and the the last term in eq.(19) r}s{l{Va)) vanishes). We will com- 
pute the self energy at finite temperature 5T,{Vg,1{Vc))- For the non-interacting ground 
state \F > with the electronic density ne(VG) we have at a temperature T the equation: 
rieiVc) = - I ^ZZ-nHVa) fF.D r^'VsT^^^ ] ' The effect of the interactions will 
replace the ground state |F > by the renormalized ground state |G >. The ground state 
represents a shifted Fermi Surface given by the self energy 5T,{Vg, 1{Vg)) determined by the 
self consistent equation: 

6E(Va,l(Va)) = 2,MVa)) < G|4|G 2,„(i(y„))4 [£M±i5(lMM^MI)] 

m=—np(VG) 

(36) 

The solution for 5E(Vg, ^(Vg)) is obtained once we replace the sum by an energy integration 
(the density of states cancel the velocity): 

6EiVa, liVa)) = hv,[{^^^^) + 7(-)F(^, ^)] • J^^^^f^fj;^ (37) 
where the explicit form e(VG, T, L) represents the effective dielectric function given by: 

i(Va,l,T,L) = 1 + - gW.)) ^ J^) . r(T)] (38) 

where r(T) represents the thermal correction, which is 1 when we use the approximation : 
IeHVG)-SvGMVG)) defF.D. r~K°?^ ] ~ 0- When the self energy is small with respect the Fermi 
energy ^^^^^(^yJ^^^ « 1, we expand the Fermi Dirac function with respect (5E(Vg, ^(Vg)) 
and find: r(T) = /^.^.[^£HI|h_fMom] _ /^^^ [£HM_«)])]. 



[1] K.J.Thomas, J.T. Nichols, N.J.AppIeyard, M.Y. Simmons,M.Pepper,D.R.Mace,and 
D.A.Ritchie Phys.Rev.Lett. 77,135 (1996) 

24 



[2] K.J.Thomas, J.T. Nichols, N.J.Appleyard, M.Y. Simmons,M.Pepper,D.R.Mace,W.R.Tribe,and 

D.A.Ritchie Phys. Rev.B. 584826(1998) 
[3] R.de. Picciotto et al. Phys.Rev.B 72,033319(2005) 
[4] S.M. Cronennwctt ct al.Phys.Rcv.Lctt. 88,226805(2002) 

[5] Y. Yoon, L. Mourokh, T. Morimoto, N. Aoki, Y. Ochiai, J. L. Reno, and J. P.Bird 

Phys.Rev.Lett 99,136805 (2007) 
[6] Y.Komijani et al. Cond-Mat/0908.2360(2009) 

[7] H. Bruus V.V. Cheianov and K. Flensberg . Physica E 10,97 (2001). 

[8] D.J. Reilly Phys.Rev.B 72,033309 (2005) 

[9] K.Hirose et al. Phys.Rev.Lett. 90,026804(2004) 
[10] D.Schmeltzer et al. Philos.Mag.B 77,1189(1998) 
[11] O.P.Sushkov Phys.Rev.B 67,195318(2003) 
[12] K.A. Matveev Phys.Rev.B 70,245319(2004) 
[13] K.Aryanpour and J.E. Han Phys. Rev. Lett. 102, 056805 (2009) 

[14] J.S. Meyer and K.A. Matveev Cond-mat 0808.2076 and Journal of Cond. Matt. 

21,023203(2009). 
[15] G.A.Fiete Reviews of Modern Physics ,79 801(2007) 
[16] V.V. Cheianov and M.B.Zvonarcv Phys.Rev.Lett. 92 176401 (2004) 
[17] O.F. Syljuasen Phys.Rev.Lett.98 textbf98,166401 (2007) 
[18] T.Rejec and Y.Meir,Nature 442,900(2006) 

[19] F.D.M. Haldane,J.Phys.C 14,2585(1981) and Phys.Rev.Lett. 47,1840(1981) 
[20] C.L. Kane and M.P.A. Fisher Phys.Rev.B 46,15233(1992) 

[21] I.Safi and H.J.Shulz Phys.Rev.B. 52, R17040 (1995) 

[22] D.Schmeltzer ,Phys.Rev.B. 63 , 1253321 (2001) and Phys.Rev.Lett. 85,4132 (2000). 
[23] D.L. Maslov and M.Stone Phys.Rev.B 52, R5539(1995) 

[24] M.Malard, D.Schmeltzer and A.Kuklov to be published in the Proceedings of the '"Interna- 
tional Conference on Strongly Correlated Electron System'" SCES(2008) 
[25] H.J. Shultz Phys.Rev.Lett. 71,1864 (1993) 

[26] D.Schmeltzer et al. Phys.Rev.B 71,045429(2005) and D.Schmeltzer et al. Phys.Rev.Lett 
95,06880(2005) 

[27] Philipe Di Francesco, Pierre Mathieu and David Senechal '"Conformal Field Theory'" page 

25 



390, Springer Verlag(1996) 

[28] N.Kawakami and S.K. Yang J.Phys: Condens.Matter 3,5983(1991) 

[29] D.Schmeltzer Phys.Rev.B 43,8650, 1991;D.Schmeltzer and A.R.Bishop Phys.Rev.B 45 
3168,1992. 

[30] D.Schmeltzer, cond-mat 0211449,vl(2002) 

[31] P.Sun and D.Schmeltzer,Phys.Rev.B ,61,349(2000) 

[32] R.Shankar ,Int.J.Mod.Phys.B 4,2371(1990) 

[33] J.M.Kosterlitz and D.Thouless, J.Phys. C5,L124(1972) 

[34] V.L.Berezmskii,Sov.Phys.JETP 34,610(1972) 

[35] A.Luther and V.J. Emery ,Phys.Rev.Lett.33,589(1974) 

[36] J.Phys.Soc.Jpn.65,3604(1996);M.Mori,M.O.Gata and H.Fukuyama cond-Mat/970802(1997) 

[37] D.Schmeltzer in preparation. 

[38] M.Pepper and J.Bird J.Phys:Condens. Matter 20,160301 (2008) 

[39] J.Rech and K.A. Matveev J.Phys:Condens. Matter 20,164211 (2008) 

[40] Y.Meir J.Phys:Condens. Matter 20,164208 (2008) 



26 



